Chiral currents in gold nanotubes 
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Results are presented for the electron current in gold chiral nanotubes (AuNTs). Starting from 
the band structure of (4, 3) and (5, 3) AuNTs, we find that the magnitude of the chiral currents are 
greater than those found in carbon nanotubes. We also calculate the associated magnetic flux inside 
the tubes and find this to be higher than the case of carbon nanotubes. Although (4,3) and (5,3) 
AuNTs carry transverse momenta of similar magnitudes, the low-bias magnetic flux carried by the 
former is far greater than that carried by the latter. This arises because the low-bias longitudinal 
current carried by a (4,3) AuNT is significantly smaller than that of a (5,3) AuNT. 
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Nanotubes and nanowires are of interest, not only be- 
cause of their potential for deployment as interconnects, 
p-n junctions and rectifiers [l], 0| in future nanoscale cir- 
cuits, but also because they exhibit fundamental physical 
properties, such as conductance quantisationy, 14 , and 
magic numbers reflecting structural stabilities j5, 6]. One 
ubiquitous property associated with nanotubes is chiral- 
ity, which arises because there is an infinite number of 
ways of rolling up a two-dimensional periodic lattice to 
form a cylinder. In addition to widely-studied carbon 
nanotubes 0], chiral nanotubes have been formed from 
a ran ge o f other materials, including gold H, B, B , plat- 
inum 10], silver 



nitride 



alkaline metals [12Ml4j and boron 
These experimental observations have 
been supported by a range of theoretical investigations 

IB-ill. 

It has recently been noted that the presence of intrin- 
sic chiral electron currents in chiral nanotubes can be ex- 
ploited to yield photogalvanic effects in heteropolar nan- 
otubes [1711. a new drive mechanism in carbon-nanotube 
windmills [23] and to produce internal magnetic fields in 
carbon nanotube solenoids [24| . In each of these exam- 
ples, the underlying lattice is hexagonal, with two atoms 
per unit cell. In contrast nanotubes formed from gold, 
silver and platinum are derived from triangular lattices 
with one atom per unit cell and therefore it is of inter- 
est examine whether or not these effects are enhanced or 
diminished compared with their carbon counterparts. 

In this paper, to answer these questions, we examine 
chiral currents in gold nanotubes. Our choice of gold is 
in part motivated by the fact that chiral currents are ex- 
pected to scale with the Fermi velocity of the underlying 
two-dimensional lattice, which in gold is approximately 
double that of graphene. 

A nanotube formed from a 2D lattice with periodic 
boundary conditions can be described by a chiral vector 
C = nai + m&2 , which defines the circumference of the 
nanotube, where ai, a 2 are the lattice vectors and n,m 
are integers. The axis of the nanotube lies parallel to the 
longitudinal translation vector T (which is perpendicular 
the the chiral vector), whose magnitude is equal to the 
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q = N , N + 1, N + 2, . . . , N + N - 1 



length of the nanotube unit cell, along the tube axis. To 
understand the currents carried by such a nanotube, we 
first calculate the electron group velocity components in 
the chiral and longitudinal directions. These are given by 
Hvc = dE(k)/dkc and Hvt = dE(k)/dkT, where kc — 
k.C and hp = k.T are components of the wavevector k 
parallel to the unit vectors C and T respectively. In the 
presence of periodic boundary conditions, kc is quantized 
and takes the values 
2nq 

where N is the number of the mini-bands and C = |C|. 
Since q and q + N are equivalent, Nq is an arbitrary 
integer. In what follows Nq is chosen such that if N is 
even, jV = 1 - N/2, whereas if N is odd, A^o = 1 =^ L . We 
also choose N to equal to the number of unit cells of the 
2D lattice contained within a unit cell of the nanotube 
and therefore — ^ < k T < fr , where T = |T|. 

The simplest model describing electronic properties of 
a triangular lattice consists of a tight-binding Hamilto- 
nian with one orbital per atom. The dispersion relation 
of such an infinite sheet takes the form 

E(k) — — 27(cosk.ai + cosk.a2 + cos k. as), (1) 

where 7 is the nearest neighbour coupling, ai = 
— §(l,v3) and a 2 = a (1,0) are lattice vectors of the 
triangular lattice, a3 = ai + a2 is an auxiliary vector 
and a is the lattice constant. For a triangular lattice 
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m— 2n 



T — ai + ^r^a 2 , where d = GCD{2m-n, 2n-m). 
For the case of a nanotube formed by imposing peri- 
odic boundary conditions, E(\s) = E(k^ ,hr) and N — 



2(n 2 



i)/d. To go beyond this simple model, 



we performed DFT calculations [26|, [27[ on AuNTs and a 
triangular gold lattice, whose lattice constant (obtained 
by relaxing the size of the unit cell) was found to be 
a = 2. 73 A. By fitting the resulting band structures to a 
third-nearest neighbour model, we obtained a more ac- 
curate dispersion relation of the form 



E(k^,k T ) 



2^7i cos(aik^ 



(9) 
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TABLE I: Fitted coupling parameters for a 2D triangular- 
lattice sheet with lattice constant a = 2.73A and for perfect 
(4, 3) and (5, 3) AuNTs (in eV). 



spect to kc and kr, the longitudinal group velocity (par- 
allel to T) is found to be 



FIG. 1: (4, 3) A AuNT and a triangular sheet showing a (4, 3) 
chiral vector. Periodic boundary conditions make the two 
lattice points at the ends of C vector identical. The T vector 
shows the unit cell of the AuNT. The red lines shows the 
1st, 2nd and 3rd nearest neighbours of the highlighted lattice 
point. 



where on = C.ai, /3j = T.ai and a4 = ai + a3, as = 
a 2 + a 3 , a 6 = a 2 - ai, a 7 = 2a 1: a 8 = 2a 2 , a 9 = 2a 3 
auxiliary vectors. 

The values obtained for these couplings are shown in 
Table HI To account for curvature effects, nine different 
couplings are needed for AuNTs, whereas, due to sym- 
metry, only three are needed for a flat sheet. With this 
choice of parameters, the ordering of the bands in eq (2) 
follows that obtained from DFT. In contrast, we found it 
impossible to obtain the correct ordering using the simple 
nearest-neighbour model of eq. (1). 

In what follows, we focus on the (5,3) and (4,3) 
AuNTs, since these are realistic experimental targets 
As a consequence of curvature, we also expect 
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that the Fermi energy of the tube will be shifted from 
that of the sheet. This shift is taken into account by an 
appropriate choice of eo- Furthermore, since the effect of 
curvature depends on the choice of n, m, the couplings 
are also allowed to vary with n, m. The corresponding 
band structures are shown in Fig. [2j 





(a)AuNT (4, 3) 



(b)AuNT (5, 3) 



FIG. 2: The black dashed lines show the DFT calculated 
band structure of the perfect folded AuNT. The red curves 
show the fitted TB band structure. Ef is the Fermi energy. 



(3) 



and in the transverse group velocity (parallel to C) is 



hv^ (fc T ) = 2 7i ai sm{ ai k^> + fck T ) (4) 
i=i 

As a reference velocity we note that the Fermi velocity 
for the triangular sheet(which varies by ±10% around the 
Fermi surface) has an average value of vf = 1.8 x 10 6 to/s, 
when averaged over the Fermi surface. 

For an infinitely-long AuNT, we now compute the 
chiral velocities of right-moving electrons, (i.e. with 
Vjf (kr) > 0), by first inverting eq. ([2]) to obtain kx as 
a function of E and q. We denote this inverse k^ + (E), 
where + sign refers to solutions belonging to branches 
with > 0. Real values of this function arise in the 
energy range £< < E < e> , where is the bottom 
of positive- VjP branch of the qth mini-band and e> is 
the top of the positive-slope branch of the gth mini-band. 
Substituting k^ + (E) into eq. ^ yields the chiral veloc- 
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ity v^ + (E) 



v, 



(9) f !.(?) + 



(kj V (E)) belonging to right- moving 
electrons of energy E. The total chiral velocity for the all 
right-moving electrons of energy E is obtained by sum- 
ming up the chiral velocities for each mini-band with a 
real k ( ^ )+ {E), to yield 



,(tot)H 



(E) 



No+N-l 



(E)0{E 



^)e{e^-E) (5) 



q=N 



To invert eq. ©, we note that since fiikx < 2tt^, 
and for most of the n, m pairs, T is much longer than a, 
the cosine functions can be approximated by parabolae, 
from which one can easily obtain an expression for kx as 
a function of energy, along with the allowed miniband en- 
ergy ranges. To obtain an explicit form for the energy de- 
pendence of the transverse velocity near an energy band 
minimum, we note that in the sum on the right-hand-side 



By differentiating the dispersion relation @ with re- of eq. if the condition < E < is satisfied for 
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FIG. 3: The top, second and third rows show chiral and lon- 
gitudinal velocities and their ratios for each miniband of the 
(4,3) and (5,3) AuNTs. The red curves show those chiral ve- 
locity values belonging to right moving electrons. Those of 
left-moving electrons are shown by dotted curves. The thick 
black curve shows the sum of the corresponding red curves, 
and the chiral velocities (red lines) for each channel. The 
velocities are in units of average Fermi velocity vf of the tri- 
angular lattice sheet. The capital letters label individual open 
channels. 



a given q, then usually it is not satisfied for —q, because a 
mini-band with no local extrenium in the Brillouin zone 
satisfies e^ 9 "* = e^ F9 ' > . However if it has local extremum 
in Brillouin zone, then both q and — q give contributions 
to the sum because if the miniband has a minimum, then 



.(-<?) 



and if it has maximum then e 



(<i) 



-(-<?) 



Therefore in the case of a minimum, the beginning of 
the energy ranges of the minibands overlap and in case 
of maximum the end of the energy ranges overlap. Con- 
sequently channels open (close) in pairs when there is a 
band with local minimum (maximum). Just as a channel 
pair opens (closes) they give a combined contribution to 
the right- hand-side of eq. (JSJ) - Noting that faKr <C 2tt 
and Taylor expanding ((3J) in kx around zero, one can 
find that the contribution from the channel pair q, —q 
is a square-root function of E. The q = case, which 
belongs to the first channel and does not have partner, 
yields a contribution 
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FIG. 4: The forces against bias voltage. The red curve belongs 
to (4, 3) AuNT and the blue curve belongs to the (5, 3) AuNT. 
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FIG. 5: The magnetic field against bias voltage. The red 
curves belongs to (4, 3) AuNT and the blue curve belongs to 
the (5, 3) AuNT. 

The above square root behavior near the bottom of a 
miniband is also found in CNTs [2J|- For gold (4,3) and 
(5,3) AuNTs, the top row of Fig. [3] shows chiral veloc- 
ities of individual minibands (red curves A-D) as func- 
tion of energy, along with the total chiral velocity (black 
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curves), obtained by summing the individual velocities. 
The second row of Fig. [3] shows the longitudinal veloc- 
ities of individual minibands. The sign of the chiral ve- 
locities alternates with successive open channels, leading 
to oscillations in the total chiral velocity with energy. To 
illustrate the behaviour of these quantities when bands 
open and close, results are plotted over a wider energy 
range than that used in Fig. [5] However results are only 
meaningful at low energies and therefore when predicting 
forces and magnetic fields in figures Fig. [4] and Fig. [5j 
we revert to energies within 1 eV of Ep. 

The chirality of right-moving electrons can play an im- 
portant role in driving nano motors, because they can 
exert a torque on the AuNT. An estimate of the maxi- 
mum possible tangential force is given by the total flux of 
tangential momentum associated with right-moving elec- 
tron injected into a bias- voltage window — ¥ > This 
takes the form 



2m 



F c = ^ I v ( c° t)+ (E)dE 

_ N a +N-1 + ~" 
q=N 



x Q(E - e!? )+ )e(4 9)+ - E)dE 

Fig. |4] shows the resulting force for (4,3) and (5,3) 
AuNTs. As expected, this is higher than in CNTs, be- 
cause chiral velocity at the Fermi energy is finite and 
also, because the Fermi velocity of AuNTs (which sets 
the velocity scale) is greater than that of CNTs. 

A further consequence of the chirality is the presence 
of an induced magnetic field, given by [24 1 
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The integrand of this expression involves the ratio of 
the chiral to longitudinal velocities. For individual mini- 
bands, these ratios are shown in the third row of Fig. 
[3l The resulting magnetic fields are shown in Fig. [5] 
Comparison between Figs. [4] and [5] shows that although 
(4,3) and (5,3) AuNTs carry transverse momenta of sim- 
ilar magnitudes, the low-bias magnetic flux carried by a 
(4,3) AuNT is far greater than that carried by the (5,3) 
AuNT. This arises because the velocity ratio of the for- 
mer is significantly higher than the latter, even though 
they possess similar chiral velocities near the Fermi en- 
ergy. 

Wc have calculated the chiral velocities carried by elec- 
trons in infinitely-long (4,3) and (5,3) AuNTs, which as 
shown in ref. [24|, can be a guide to the size of chiral 
currents in finite NTs connected to reservoirs. We have 
found that in a similar fashion to CNTs, chiral currents 



are oscillatory functions of energy, but unlike in CNTs, 
the chiral current has a finite value near Ep. Further- 
more the tangential force and the induced magnetic field 
is higher than in comparable CNTs. We have consid- 
ered perfectly-periodic nanotubes only. For the future it 
will be of interest to consider the effect of disorder on 
chiral currents. In this regard, one notes that at least 
in one dimension, disorder, which preserves the average 
spatial symmetry of a lattice, does not completely ran- 
domize the phase and density of current-carrying states 
28l 2§|. Therefore, one expects chiral currents to persist 
in the presence of disorder, although with a diminished 
magnitude. 
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